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1 INTRODUCTION 

Quantum chromodynamics (QCD), the theory on strong interactions, has been 
tested extensively in hard, large-momentum processes. In these processes, a large 
amount of energy or three- momentum is transferred to one (e.g., deep inelastic 
scattering) or few (jet production) quark and gluon constituents of initial hadrons. 
Owing to the asymptotic freedom property of QCD, the coupling strength is 
small. The time scale for this part is a small fraction of the time scale of the 
overall process. These two properties allow the factorization and perturbative 
treatment of the hard part of the process from the rest of the matrix element. 
This programme has been very successful, and no clear discrepancies have been 
found between experimental results and theoretical calculations. 

The goal of the experimental heavy ion program at ultra-relativistic energies 
is to study QCD in an environment very different from that encountered in hard 
processes, in a dense system of quarks and gluons. When the heavy ion pro- 
gramme started two decades ago, the original goals were the production and 
study of dense, thermally equilibrated, strongly interacting matter, the quark- 
gluon plasma (QGP). Although this is still the highest priority, phenomena other 
than the formation of thermally equilibrated QGP also occur in dense partonic 
systems and can be studied in heavy ion collisions. E.g. the existence of what 
is called a Color Glass Condensate (CGC), which is related to the saturation of 
gluon occupation numbers in dense components of initial wave functions, may 
be important for the formation of the QGP and may also have other observable 
effects (1). Also, the phenomenon of jet quenching, the loss of energy of a high- 
energy parton (quark or gluon) when it traverses a high-density parton system, 
can be expected to occur in the dense environment of heavy ion collisions, even 
when the system is not fully equilibrated. 

The earlier experiments at Brookhaven National Laboratory (BNL) and CERN- 
SPS provided clear evidence of collective phenomena in nuclear collisions. Re- 
cent results from experiments at the Relativistic Heavy Ion Collider (RHIC) 
at BNL show that total multiplicities exceeding 1000 particles per unit rapid- 
ity, dN/dy > 1000, are produced in central (head-on) collisions of gold nuclei 
at center-of-mass energies up to ^snn = 200 GeV. Here ^Jsnn is the center- 
of-mass energy for a nucleon pair. The measured average transverse momen- 
tum is pt ^ 0.5 GeV, indicating a total energy per unit rapidity interval of 
dEx/dy > 500 GeV. At time r after the nuclei have passed through one an- 
other, the volume occupied by the produced quanta in a rapidity interval Ay is 
AV = rAyAx, or tAt = tttR 2 a for unit rapidity. At r = 1 fm/c the above 
numbers imply (2) a particle density of ~ 10 fm -3 and an energy density of 
~ 5 GeVfm -3 , well above different estimates and the results from lattice calcula- 
tions of the energy density at which the phase transition from confined hadrons 
to unconfined quarks and gluons occur. 

The subject of this review is the use of hydrodynamic modeling to describe 
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the expansion and dilution of matter produced in nuclear collisions. One motiva- 
tion for this is to formulate a framework to study different observable quantities 
and correlations among them. Obviously hydrodynamics alone does not suffice 
because, first, at high energies, particle production can not be included and, sec- 
ond, only the properties of produced particles, not the hydrodynamic densities 
during the expansion, can be measured directly. Below we describe a possible 
dynamic scenario to calculate the production of initial matter. This provides the 
initial conditions for solving the hydrodynamic equations. Because the produc- 
tion dynamics is still not completely under control, sometimes it may be useful 
to use physically motivated parametrizations of initial conditions to study, for 
example, how the details of particle spectra depend on different features of initial 
conditions. The hydrodynamic description also needs another supplement. A 
link from hydro dynamical quantities to particle spectra is necessary at the end 
of expansion, when the particles become independent and fly to detectors. 

The use of hydrodynamic concepts like temperature, pressure, and flow velocity 
cannot be strictly justified for matter formed in a heavy ion collision. Although 
the total number of particles produced is several thousand at RHIC and may be an 
order of magnitude more in the future ALICE experiments at the Large Hadron 
Collider (LHC) at CERN, they hardly form a macroscopic system in the proper 
thermodynamic sense. However, one can argue for partial equilibration and the 
formation of collective phenomena from the numbers given above. At initial 
particle densities of the order of 5-10 fm -3 , and even with modest estimates of 
cross sections of the order of 1-2 mb, mean free paths are < 1 fm, much below the 
nuclear size 2Ra > 10 fm. Therefore, frequent collisions occur and momentum 
is transferred from denser regions toward less dense regions. In describing the 
main consequences of these secondary collisions, concepts like temperature and 
pressure are useful. 

Although the use of hydrodynamics may be justified for some features of nuclear 
collisions, that need not be the case for others. Defining energy density and 
pressure in the usual way in terms of the local momentum distribution, one 
obtains for pressure of massless particles P = e/3 for any isotropic momentum 
distribution f(p), p = \p\ = E. In this case, describing the build-up of collective 
motion or flow using hydrodynamics with P = e/3 as the equation of state (EoS) 
could be a reasonable approach, even when the form of p dependence of f(p) 
differs from that of the thermal equilibrium distribution, fth(p)- However, if f(p) 
differs significantly from /th(p), any conclusions based on the detailed momentum 
dependence of f(p) would fail. This could well be the case at large momenta. 
Initial production is expected to contain a component that has an approximate 
power-law behavior at large momenta. It will take longer and requires a larger 
volume than available to change this power behavior to the exponential form 
of thermal distribution. The observed behavior of high-px hadron spectra shows 
clearly that the high-momentum partons are not thermalized and that they suffer 
an energy loss while traversing the produced dense matter. Because the high- 
energy partons form only a small part of produced matter both in multiplicity 
and in transverse energy, a thermal equilibrium description can still be adequate 
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for the bulk of the matter. 

Keeping in mind all the reservations, we review the use of hydrodynamics in 
describing the heavy ion collisions at collider energies, and the calculation of 
observable quantities. We compare the calculated results to the measurements 
mainly at RHIC, mentioning some results from the CERN-SPS. We also present 
examples of predictions for the future ALICE experiment at CERN. 

2 HYDRODYNAMIC EXPANSION 

Hydrodynamics is the theoretical framework describing the motion of fluid, a 
continuous, flowing medium. The equation of motion can be derived from kinetic 
equations. Hydrodynamic equations take the simplest form if local thermal equi- 
librium is assumed. In this treatment, there are no dissipative effects. To take 
such effects into account by approximation, small deviations from the local equi- 
librium are assumed. A linear treatment of deviations leads to a system of equa- 
tions which contains the viscous coefficients including heat conductivity, in the 
case of conserved currents. For a review, see Reference (3). In most applications 
of hydrodynamics to heavy ion collisions, viscosity has been neglected. In studies 
with viscosity, results on global, integrated quantities do not differ qualitatively 
from those without viscosity. However, for example, transverse spectra at larger 
transverse momenta px > 1-5 GeV may start to deviate clearly from ideal-fluid 
calculation (4-6). One should remember, though, that the viscous properties of 
strongly interacting matter are not well understood, and the approximations in 
the numerical work also introduce uncertainties. 

Once EoS is known and initial conditions are specified, the hydrodynamic equa- 
tions determine the expansion of the fluid. In the context of describing heavy 
ion collisions, the use of these equations requires knowing the EoS of strongly 
interacting matter and the primary production of particles. Detailed knowledge 
of microscopic processes is not required if a very strong assumption is taken: The 
expanding system stays in local thermodynamical equilibrium. This becomes of 
great practical importance if one wants to include in the hydrodynamic expansion 
the transition from quarks and gluons to hadrons. The complicated deconfine- 
ment or hadronization processes need not be known in microscopic detail; all that 
is necessary is the thermodynamic equation of state as computed, for example, 
in lattice QCD. 

2.1 Hydrodynamic Equations 

The hydrodynamic equations 

8^ = 0, T» v = {e + P)u»u u -Pg»" (1) 

express, in terms of the energy- momentum tensor T^ v , the conservation of energy 
and momentum in continuous, flowing matter. The quantities defining T^ v are: 
e := the energy density, P := the pressure, and vJ 1 := the flow four-velocity, 
normalized to u^u^ = 1 as usual. The simple form of T^ v above holds for an 
ideal fluid. 
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In addition, if the system contains conserved densities rij, such as those of 
charge and baryon number, their evolution is expressed by continuity equations 

M = • (2) 

With several conserved densities or non-zero viscous terms, the definition of is 
not unique. One can define the velocity in terms of one of the conserved currents 
by writing = riiU^, where n-i = y/ji^jf, or one can use the energy-momentum 
tensor to define the flow velocity In the first case, usually referred to as the Eckart 
definition, there is no flow of charge Qi in the local rest frame uf = (1,0), the 
Eckart frame. For a non-ideal fluid, the energy flow would usually be non-zero 
in the Eckart frame. A definition of the fluid velocity in terms of the energy- 
momentum tensor, referred to as the Landau definition, is such that the energy 
flow is zero in the local rest frame, but usually the flow of different charges does 
not vanish. Here, the only charge we consider is the baryon density. We also 
treat the matter in the final state as an ideal fluid, and thus the two choices for 
the flow velocity coincide. For a review, see Reference (3). 

The only properties of the dynamics contained in Equations 1 and 2 are the 
conservation laws. However, the relations between the thermodynamic densities 
e, P and n«, or alternatively their definitions in terms of temperature T and 
chemical potentials Hi — e.g., e = e(T, //j) — constituting the EoS, depend on 
the details of the dynamics among the constituents of the matter. The need for 
an EoS is obvious: The Equations 1 and 2 contain five equations whereas there 
are six quantities to be defined by solving the equations. These quantities are the 
three components of the velocity, the energy density, pressure and the baryon- 
number density. For a non-ideal flow, transport coefficients would enter into the 
expressions of the energy-momentum tensor and currents, and their derivation 
from theory requires the knowledge of microscopic dynamics in the same way as 
the derivation of EoS. 

Most of the detailed hydro dynamical discussion below is limited to the situation 
of scaling longitudinal flow and invariance under longitudinal Lorentz boosts. 
This means that the longitudinal flow velocity is v z = z/t, and hence, the flow 
rapidity is rf = log[(i + z)/(t — z)], which is also often termed the space-time 
rapidity. If the initial densities are assumed to depend on t and z only through 
the longitudinal proper time r = y/t 2 — z 2 , e.g., e = e(r, r), the expansion will 
evolve so that densities remain independent of rj and the v z will retain the scaling 
form v z = z/t. In this situation the system is said to be boost-invariant (2). 

Before discussing the EoS in more detail, we note that the most useful form 
of the EoS for solving the hydrodynamical equations is provided by the relations 
among the densities e, P, and rij when the hydrodynamical equations are written 
in the form of Equations 1 and 2. In this form, temperature and chemical poten- 
tials do not appear in the equations. For the calculation of observable quantities, 
such as hadron spectra or electromagnetic emission, relations among the above 
densities and the temperature and chemical potentials must be specified, as we 
discuss below. 

The other ingredient that must be provided from outside into the hydrodynamic 
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description are the initial conditions, e.g., in terms of initial energy distribution 
and velocity. From the physics point of view this is a very profound problem be- 
cause it involves not only the primary production dynamics, but also the question 
of thermalization of produced particles. We discuss first the initial conditions, 
then the EoS, and finally the calculation of physical observables. 

2.2 Initial Conditions 

As mentioned above, primary particle production cannot be formulated within 
the hydrodynamic framework in a realistic way in high-energy nuclear collisions. 
The dynamics of particle production is a separate problem and, if solved, it pro- 
vides the initial conditions for the hydrodynamic expansion. Initial conditions 
specify the thermodynamic state of the matter and its velocity on an appropri- 
ate space-time boundary, which, for example, in the boost-invariant case dis- 
cussed below, can be taken to be a constant, longitudinal proper time surface 
r = \Jt 2 — z 2 = To- From the point of view of the hydro dynamical calculation, 
the initial conditions can be provided either by a dynamical calculation of pri- 
mary particle production or by a reasonable parametrization, with the parameters 
either given by physical arguments or fixed by comparing (some of) the results 
with experimental data. 

There are different approaches to primary production, such as pQCD + final- 
state saturation (minijet) (7) and the color glass condensate model (8) based on 
the initial state parton saturation (9). Both describe the produced matter as a 
parton system. Also, models based on string formation and decay, such as the 
DPMJET model (10), are used for the calculation of final hadron spectra. In such 
a model, a varying fraction of energy is in the form of color strings at an early 
stage of the collision, and they are not readily connected with a hydrodynamic 
description that assumes particle- like constituents of matter. We do not consider 
string models further, but we describe briefly parton-based approaches to the 
primary production. 

To illustrate different key factors that enter the determination of initial condi- 
tions from dynamical calculation of particle production, we consider a perturba- 
tive QCD calculation of parton production as an example. For such a calculation 
to converge, a cut-off must be provided on small momentum-transfer collisions. 
In this model the cut-off is obtained from a saturation condition expressed in 
terms of the transverse nuclear geometry and the number of produced partons. 
At collider energies, the saturation scale turns out to be typically p S at ~ 1 ... 3 
GeV, which is clearly larger than Aqcd ~ 0.2 GeV. Because this cut-off is smaller 
than what is usually used in jet calculations, produced partons are often called 
minijets, as partons close to the cut-off dominate the production. In addition to 
the cut-off, the ingredients of the calculation are the parton distribution functions 
of colliding nuclei and the parton-parton cross sections. These cross sections can 
be calculated from basic QCD theory, but parton distributions must be provided 
from other measurements. 

The nuclear parton distributions are usually expressed in terms of parton dis- 
tributions of nucleons that, however, are known to be modified in nuclei. The 
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nuclear modification factor R A (x,Q 2 ) is the nuclear parton distribution normal- 
ized to a single nucleon and divided by the parton distribution of a free nucleon. 1 

The perturbative QCD calculation of minijet production is a momentum-space 
calculation, as is the case in most production models. To define the initial spatial 
densities, a connection between the momentum of a produced parton and its 
space-time formation point is needed. At collider energies, the hard partons of 
the colliding nuclei are Lorentz contracted to a region on order of 2R A / 7cm < < 1 
fm. We consider the collision region as a point in the longitudinal direction that 
allows us to assume that the rapidity of the minijet coincides with the space-time 
rapidity of the formation point, y = rj = (l/2)ln[(i + z)/(t — z)\. We take the 
formation (proper) time to be the inverse of the saturation scale, tq = l/p sat . 
Thus, the minijet matter forms along the hyperbola t = y 7 z 2 + Tq with initial 
longitudinal flow velocity v z {tq) = z/t. To determine the transverse distribution, 
we must start with a calculation of production cross section. 

To obtain the initial conditions for baryon-number density and energy density, 
we first need the minijet cross sections for (anti)quarks and gluons and their first 
moments in transverse energy (momentum) in nucleon-nucleon collision, each 
calculated in a rapidity interval Ay and integrated in px from the saturation 
cut-off pt = Psat to its maximum value (11, 12): 



Cjetfeat, Vs, Ay, A) = / dp T ^ , 

IT?\( rA A\ f^A da iet (^,Ay,A) 
Ojet {E T ) (p sa t , Vs, Ay, A) = dpT Pt — -, 

J Psat ^PT 



'Psat 

The total number of minijets and the total amount of transverse energy in Ay in 
a nucleus-nucleus collision is obtained by multiplying the corresponding nucleon- 
nucleon cross section with the nucleon-nucleon luminosity of the collision (in- 
cluding an extra factor of two for the number of minijets). This is given by the 
overlap function T^s(6) of transverse densities Ta(b)( s ) of the colliding nuclei: 

T AB (b) = J d 2 sT A (\b - s\)T B (s) = T AB (b) , 

/+oo 
dz PA (z,s) =T A (s), 
-oo 

where b is the impact parameter and s the transverse coordinate in nucleus A. 
For example the number of partons (which can be defined only in lowest order 
because it is not an infrared-safe quantity at higher orders) produced in a central 
zero-impact-parameter collision of equal nuclei is 

AN AA = T AA (0)aj ct (p sat , Va, Ay, A) (3) 

in a rapidity interval Ay. A similar expression with <7j ct replaced by aj e t(-Er) 
gives AEt, the transverse energy of minijets, in Ay. 

Before discussing how to formulate the saturation condition to fix p sat , we no- 
tice that the average densities are obtained by dividing the total quantity with the 

1 In the actual calculation protons and neutrons are treated separately. 
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volume that corresponds to the rapidity interval Ay, AV = AzAt = T$AyT:R 2 A . 
This procedure, with densities averaged over the transverse plane, is easily gen- 
eralized to local densities. The nucleon-nucleon luminosity in a transverse-area 
element d 2 s is Tk(|b — s|)Tb(s)oPs, and the volume element is dV = Azd 2 s = 
rAyd 2 s, leading to (13) 

"^ D(T - s) = ^s = 5; r - <(|b - s|>rB(s) (4) 

for the parton density and 

£ ^ s > = ^ = W IM|b - s|,r * (s) < 5 > 

for the energy density. The densities depend on the cut-off scale through the 
cut-off dependence of cross sections. 

The minijet cross sections above can be calculated separately for gluon, quark, 
and anti-quark jets, allowing for the separate determination of the densities of 
quarks and anti-quarks. From these densities the initial net baryon number den- 
sity is obtained as = (n g — rig)/ 3, which provides the initial condition for the 
net-baryon- number current that satisfies the conservation law (Equation 2). 

Up to here, we have essentially discussed how to obtain from a boost-invariant 
momentum-space calculation of production cross sections, <7j ct and <7jet(-Er), the 
local densities n# and e. To close the calculation of minijet cross sections, the 
saturation momentum p sat must be specified, and we do this by assuming that the 
parton (mainly gluon) production saturates when the wave functions of produced 
partons start to overlap. 

In the transverse direction the scale of the wave functions is l/p sa ,t- The scale 
in the longitudinal direction is not as obvious, but we assume it is the same at 
the production time To = l/p sa .t- At this time, the volume occupied by particles 
in the rapidity interval Ay equals AV = ToAyAj-- At saturation, dividing this 
volume with the volume occupied by one jet, Vj e t, should equal the number of 
produced partons ANaa in the rapidity interval Ay. This leads to the condition 

AN A A(Psat,y/s,Ay,A) vr 2 



Ay (p. 



sat) 



\2 



= A T = ttR a . (6) 



To avoid introducing a rather arbitrary cut-off parametrization at the nuclear 
edges in transverse plane, we do not try to define the calculation of p sat locally in 
the transverse plane. Instead, an effective value p sa ,t, obtained from Equation 6 
above, is used. Solving p sa t from this equation completes the calculation of pri- 
mary production in the pQCD+saturation model. The initial energy densities at 
l/fsat at RHIC and LHC energies, ^/J/vtv = 200 and 5500 GeV, respectively, are 
shown in Figure 1. 

The above formulation of calculating the initial densities is applied in a central 
rapidity bin \y\ < Ay. In the pQCD calculation, particle production depends 
on the rapidity through the parton distributions of colliding nuclei. However, 
when we use the results as initial conditions for the hydrodynamic calculation, 
we assume that boost-invariance is a good approximation at y pt and take the 
densities to give the boost-invariant initial conditions with no rj dependence. 




Figure 1: Transverse dependence of the initial energy distribution for a gold-on-gold 
collision at the Relativistic Heavy Ion Collider (RHIC) (dashed line) and lead-on-lead 
collision at the Large Hadron Collider (LHC) energy (solid line). The saturation scale 
is ftat = 1-16 GeV at RHIC and 2.03 GeV at the LHC, with formation times of 0.170 
fm/c and 0.100 fm/c, respectively. The dashed-dotted line shows the energy density at 
the LHC if Tj = 0.170 fm/c is used, emphasizing the strong dependence of initial energy 
density on the assumed initial time. 



An attractive approach to particle production in heavy-ion collisions at collider 
energies has been based on the assumption that the initial-state parton densities 
saturate and non-linear dynamics becomes dominant (9). With increasing colli- 
sion energy, the small-a; partons, in particular the gluons, become the dominant 
part in the production, and their density in the initial wave function of the nu- 
cleus becomes so high that gluons interact coherently and their density saturates. 
The momentum below which gluons saturate is called the saturation scale, Q s . 
It depends on the collision energy yfs and the mass number of the nuclei. From 
the point of view of the color fields, the high density or large occupation numbers 
of the field quanta with momenta < Q s can be described as the formation of a 
color glass condensate. This part dominates the primary production at large \/s 
and can be treated in terms of a classical effective field theory (1, 8). Quantum 
corrections to the classical treatment have also been considered ((14, 15); see also 
Reference (16) and references therein). 

In the effective field theory approach to gluon production in AA collisions, when 
boost-invariance is assumed, a gauge can be chosen such that the problem can 
be formulated as a dimensionally reduced 2+1-dimensional theory. A numerical 
approach (17) with lattice regularization can be applied to the reduced theory. 

The numerical calculations with lattice regularization started using SU(2) sym- 
metry (18, 19) and other simplifications, such as cylindrical nuclei, but were soon 
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formulated for SU(3) (20, 21) and realistic nuclear geometry. Also, local color 
neutrality in the transverse overlap region of the collision (21) was imposed and 
other smaller inconsistencies corrected, leading to a ratio of E T /N, which is con- 
sistent with the physical interpretation of the saturation scale (22). 

The lattice approach does not, however, give a value of the saturation scale 
itself; the overall normalization must be obtained from elsewhere. For RHIC 
phenomenology, the authors of Reference (21) suggest two sets of results, which, 
in light of the latest results of Reference (22) lead to following, qualitatively 
different descriptions of the final state: 

• In the case of a smaller scale jj, the total transverse energy E T ~ g^R 2 A u^ 
produced from the classical fields roughly equals the experimentally mea- 
sured result, whereas the number of initially produced partons (N ~ g 2 R\n 2 ) 
is only about half of the multiplicity of hadrons measured in the experiment. 
In this case, the only change in the final state is the fragmentation of par- 
tons to ~ 2 hadrons on average. However, there would be no significant 
hydrodynamic evolution because this would reduce the transverse energy 
below the measured value. In this picture, which corresponds to the sce- 
nario suggested by Kharzeev & Levin (23), one would expect the photon 
and lepton pair emission after the primary interactions to be very rare. 

• For a larger saturation scale, the number of partons is close to the mea- 
sured number of hadrons but the initially produced transverse energy is 
approximately 2.5 times bigger than the measured one (22). In this case, 
production must be followed by strong initial collective expansion in the 
longitudinal direction, allowing for a transfer of energy into the longitudi- 
nal motion. This case corresponds to the evolution suggested by pQCD + 
saturation + hydrodynamics calculation (13). 

The energy dependence in the above models enters through the dependence of 
the saturation scales on the center-of-mass energy leading, for example, to rather 
similar growth of multiplicity from RHIC to the LHC in both models. 

For a head-on, zero-impact-parameter collisions, the produced system is cylin- 
drically symmetric, all quantities depend only on r and r, and the transverse flow 
is radial with no azimuthal dependence. We show results on hadron spectra dis- 
playing effects from radial flow. However, a good test for the applicability of the 
hydrodynamic description of heavy ion collisions is provided by nonzero-impact- 
parameter collisions without cylindrical symmetry. In our example of the calcu- 
lation of primary densities the expressions (4, 5) hold also for nonzero-impact- 
parameter collisions. However, the determination of saturation scale becomes 
more involved, and a straightforward generalization of the saturation condition 
(Equation 6) leads the multiplicity to have a too-flat dependence on the number 
of participant nucleons, a possible measure of non-centrality of the collision (24). 

Above, the transverse dependence of the initial densities is given by the num- 
ber of collisions per unit transverse area: n co ii(b, s) oc T^(|b — s\)Tb(s). In 
another popular phenomenological approach, the initial densities are assumed to 
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be proportional to the number of participants, also known as wounded nucleons 2 
per unit transverse area. In the eikonal Glauber model, this is defined as 



where a is the nucleon-nucleon cross section. For a zero-impact-parameter col- 
lision, this indicates, except at the edges of the nuclei, a radial dependence ap- 
proximately proportional to the sum of thickness functions Ta{t) and Tg(r). In 
the central region of the overlap area, the resulting density distribution is flatter 
than when the density is proportional to the number of collisions. This means 
smaller pressure gradients and slower evolution of transverse flow. 

In the literature, both the proportionality to the number of binary collisions 
and to the number of participants has been used to fix either the initial entropy 
density or energy density. The proportionality constant is chosen to reproduce 
the measured final particle multiplicity in most central collisions. The centrality 
dependence of the multiplicity is then predicted by the model, but at RHIC, 
neither binary-collision scaling nor wounded-nucleon scaling reproduces the data. 
However, a linear combination of them does and is therefore used to describe the 
initial density distribution (25, 26). 

We have not yet specified the initial transverse velocity vy. Usually this is 
taken to be zero. This choice is supported by the argument that the final state 
in each primary collision is randomly oriented in the transverse plane and thus 
one expects the transverse-momentum average in any volume element to vanish. 
There is a slight flaw in this argument because the transverse density of produced 
partons is not constant, and this can lead to a nonzero momentum average in 
a (finite) volume element. However, comparison to experimental data shows 
that for agreement, only small initial transverse velocities are allowed. In the 
calculations shown below, vt{to,t) = has been used. 

Because the use of hydrodynamics presumes thermal equilibrium, the time scale 
for thermalization after the primary production must be fixed. The dynamics of 
thermalization can be even more difficult to solve than that of primary produc- 
tion. Results of theoretical studies of thermalization have not yet converged. In 
the so-called bottom- up thermalization scenario (27), thermalization times are 
predicted to be long, of the order of 2. . . 4 fm/c. A more recent idea of the role of 
instabilities in the expansion predicts much shorter thermalization times of order 
below 1 fm/c (28). 

The thermalization time scale is an important issue because a hydrodynamic 
description of the elliptic flow can be achieved only if the thermalization time is 
short. In the numerical examples we provide below, we always use thermalization 
times below 1 fm/c. When showing results based on the minijet initial state, we 
use the production time scale as the thermalization time, = l/p S at- 

2 Strictly speaking wounded nucleons are nucleons that scatter inelastically whereas partici- 
pants are nucleons that scatter elastically or inelastically. In the recent literature this difference 
is usually ignored and both terms are used in the sense of participants. 
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2.3 The Equation of State of Strongly Interacting Matter 

A major complication in the description of the evolution of the matter produced 
in a high-energy collision of nuclei is the change in the degrees of freedom. The 
dense initial-state parton matter expands and turns into a gas of hadrons and 
hadron resonances when dilute and cool enough. According to the present under- 
standing, at large \xb and small or moderate T, there exist different correlated 
phases, such as the phase with color-flavor locking (29). Nearing smaller ub, the 
transition between hadron resonance gas and QGP is believed to be of first order 
when ub is not too small. For two light quarks and one heavy quark, the phase 
boundary is conjectured to end at a critical point, and below that the transition 
is a rapid cross-over. The quantitative theoretical information from QCD lattice 
simulations, which have recently been extended from the T-axis to finite chem- 
ical potential, supports this picture. There are also arguments in favor of the 
existence of strong correlations in the quark matter close to the phase boundary. 
These may explain the ideal-fluid behavior of QGP indicated by successful hy- 
drodynamical explanation of elliptic flow. When matter is assumed to be in the 
state of non-interacting quarks and gluons, the ideal QGP, a simple ideal gas EoS 
of massless particles (P = e/3) is often used to describe it. A more sophisticated 
but less usual way is to use parametrized lattice QCD results. 

In heavy ion collisions at collider energies the net-baryon number is small, with 
U B 5i 50 MeV (30), indicating a cross-over transition. Somewhat unexpectedly, 
from the point of view of hydrodynamic expansion, the difference between a weak 
first-order transition and a rapid cross-over is not very significant, as long as the 
EoSs are relatively similar away from the transition region and the increase in 
entropy and energy densities around the critical temperature is sufficiently large 
and rapid (31). Upon a closer look this is not so surprising because the main 
qualitative feature is a jump in the thermodynamical densities e and s. The 
size of the jump depends essentially on the size of the change of the number of 
degrees of freedom. From the point of view of hydrodynamics, the rapid jump in 
e combined with a slower change in pressure appears as a softening of the EoS. It 
is seen as a slow-down in the acceleration of the transverse flow in the transition 
region. Details differ for the two cases, but the final features of flow are quite 
similar and the quantitative differences in the final hadron spectra are small. 

An interacting hadron gas can be described in good approximation as a gas of 
noninteracting hadrons and resonances. The inclusion of resonances mimics the 
effects of both attractive and repulsive interactions between hadrons reasonably 
well (32). However, the repulsive interaction between baryons at large net-baryon 
densities must be included as an additional mean field (33) or as an excluded 
volume correction (34) to give a reasonable phase-transition behavior between 
hadronic and partonic phases. A detailed account of constructing an EoS with a 
mean field can be found, for example in Reference (35). 

In calculations, we use an EoS with ideal QGP in the high-temperature phase 
and a hadron resonance gas with a mean field below the transition. A first-order 
transition is implemented by introducing a bag constant B into the QGP phase 
and connecting the two phases with a Maxwell construction. We use Nf = 3, 
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and the bag constant B and mean field constant K are chosen to be i? 1 / 4 = 243 
MeV and K = 450 MeV fm 3 , respectively, giving T c = 167 MeV for the transition 
temperature. 

An additional complication in constructing an EoS of a hadron gas relevant 
for heavy ion collisions is the chemical composition of the hadron gas. The usual 
assumption of hydrodynamics is that of chemical equilibrium. This assumption is 
supported by thermal models that can reproduce the observed hadron abundances 
by assuming a thermal source in T pt 170 MeV temperature. However, many 
studies have found that the px distributions of hadrons are better described by 
assuming a colder, flowing source in T = 100-140 MeV temperature. Thus, the 
assumption of chemical equilibrium between these temperatures is questionable. 
In many hydro dynamical calculations, this observation is simply ignored and 
chemical equilibrium is assumed to hold until kinetic freeze-out at T = 100-140 
MeV. In such cases one can reproduce the slopes of the hadron px spectra, but it 
is not possible to reproduce simultaneously both baryon and antibaryon yields. 

One solution to this problem is the so-called single freeze-out model (36), in 
which a suitable choice of freeze-out surface allows one to fit the px spectra, 
even if the kinetic freeze-out temperature is taken to be the same T ~ 165 MeV 
as the chemical freeze-out temperature. As we show below, in the context of 
hydrodynamical models, similar approach with T ~ 150 MeV can be used to 
reproduce, at least approximately, the hadron px distributions in most central 
collisions at RHIC energies (37). Whether the anisotropics of particle distribu- 
tions (see section 3.3) can also be reproduced this way, has not been tested so 
far. 

Another solution to this problem is to assume two separate freeze-outs - 
chemical and kinetic — and to modify the EoS between these temperatures ac- 
cordingly. In such an approach, the hadron yields are assumed to be fixed at some 
chemical freeze-out temperature, usually soon below or at the phase-transition 
temperature. These hadron yields are subsequently described as conserved cur- 
rents, and each conserved hadron species is assigned a chemical potential. This 
way the yields of all hadron species can be reproduced separately, even if a low 
kinetic freeze-out temperature is used (38). 

Such an EoS changes the buildup of collective motion — i.e. flow — very little 
because pressure as a function of energy density, P = -P(e), is very similar to the 
chemical equilibrium EoS (39, 40). However, temperature as a function of energy 
density changes radically, and when collective and thermal motion are folded to 
calculate observable particle distributions, the results differ (40, 41). 

2.4 Transverse flow 

To illustrate the transverse flow, Figure 2 shows the boundaries of QGP and 
hadron gas, with the mixed phase between them. Also, three contours of tem- 
perature are depicted, as well as the flow lines with 1-fm intervals starting at To. 
The initial conditions are those from the pQCD + saturation model at RHIC 
energy y / ijv r /v = 200 GeV. 

According to this calculation, the maximum lifetime of the plasma phase is 
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Figure 2: Temperature contours at 300 (in quark-gluon plasma QGP), 150 and 120 MeV 
(hadron resonance gas, HRG), and the boundaries of mixed phase (MP) with QGP and 
HRG at T c — 167 MeV. Flow lines are also shown. Initial conditions are from a pQCD + 
saturation calculation at ^Snn = 200 GeV. Note that the slope of the flow line is related 
to the velocity and the curvature to the acceleration. 

5 fm/c, and that of the mixed phase is ~ 8 fm/c. The flow lines show how the fluid 
elements accelerate and move. The slope of the flow line is related to the velocity 
of the fluid and the curvature to the acceleration. In the mixed phase, where 
the pressure is constant, flow lines are straight because there is no acceleration. 
At small r, gradients are small and the flow lines bend slowly. In particular, 
in the plasma, when r increases, the pressure gradient, and consequently the 
acceleration, increases, indicated by the faster bending of flow lines. Along the 
flow line starting at r = 8 fm/c, the densities are small even at To, and this region 
is insignificant in calculating the spectra. For calculational reasons the initial 
densities are taken to go smoothly to zero, and the hydrodynamic equations are 
also solved at large values of r. 

3 HADRON DISTRIBUTIONS AND CORRELATIONS 

3.1 Freeze-out and the calculation of hadron spectra 

As matter expands, distances between particles become large, collisions cease and 
momentum distributions freeze out. The condition for the freeze-out to occur is 
usually expressed locally in terms of the energy density or temperature reaching 
a given value. This determines a three-dimensional freeze-out surface a^{x) in 
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space-time. The prescription of Cooper and Frye (42) convolutes the flow with 
the thermal motion along the freeze-out surface: 
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The second expression above is valid for cylindrically symmetric, boost-invariant 
flow with v r the radial flow velocity, "y r = — v%, and Ki and Ii are Bessel 

functions of second kind. 
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Figure 3: The effect of flow on the spectrum of kaons. Temperature is kept unchanged 
and the spectrum is shown for radially flowing matter at velocities v r — 0, 0.6 and 0.8. 



The unstable particles are treated as follows: First, the spectra of all hadrons 
and hadron resonances are calculated using Equation 9. We then follow the chains 
of all possible two- and three-body decays and collect the spectra of final stable 
hadrons (43). Stable hadrons can either be interpreted as all those that do not 
decay through strong interactions, or we can follow the feed-downs via weak and 
electromagnetic decays. E.g. we can calculate the 7r° spectrum including both 
the 7r°'s at freeze-out and all the decays that lead to 7r°'s, and then study the 
photon spectrum from 7r° decays alone or from all decays with photons in the 
final state. Phenomenologically an important case is that of feed-down nucleons 
from weak decays of strange hyperons. Here again, we can calculate, for example, 
the spectrum of A's and then see how it contributes through weak decays to the 
spectrum of protons. 
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Before showing our results at RHIC and LHC energies we discuss how the 
radial flow affects a spectrum. In Figure 3 we show the spectrum of kaons from 
matter at rest or flowing at velocities v r = 0.6 and 0.8. We assume that the 
matter decouples at a fixed time so that only the second term in Equation 9 
contributes. The temperature of the matter is the same in each case and is 
essentially the inverse of the logarithmic slope of the spectrum for v r = 0. From 
the asymptotic properties of modified Bessel functions of second kind, it is clear 
that at large wit the change in the slope, when v r changes to a nonzero value, 
can be expressed by replacing the temperature T with an effective temperature 
T e ff = Ty/(1 + v r )/(l — Vr). For v r = 0.6 the change is a factor of two and for 
v r = 0.8 a factor of three. 

3.2 Transverse momentum spectra of hadrons 

Next we compare some of the hydrodynamical results (44) with the experimental 
transverse momentum spectra measured by the STAR (45-48), PHENIX (49-54), 
PHOBOS (55, 56) and BRAHMS (57-59) collaborations for the most central bins 
in Au+Au collisions at ^snn = 130 and 200 GeV. The calculated spectra are 
obtained by using Equation 9 with the flow illustrated in Figure 2. 

Note that a hydrodynamic calculation cannot describe the hadron spectra at 
large transverse momenta. At large pr the hydrodynamical calculation shows an 
approximate exponential behavior, whereas the tails of measured spectra essen- 
tially obey a power law. At RHIC, the transition from steep exponential to a 
shallower power behavior takes place at pr ~ 3 GeV. The fraction of hadrons 
with p T > 3 GeV from all hadrons is small, and they originate from the fragmen- 
tation of high-energy partons, which suffer some energy loss in the dense medium 
of low-energy partons, but are not thermalized. We return to the interplay of the 
low-energy partons, which provide the main transverse energy and are assumed 
to thermalize, and the high-energy partons, which lose some fraction of energy 
in rescattering but require a much larger system for thermalization. 

We start with the px spectra of identified hadrons at midrapidities. Figure 4 
shows the PHENIX data collected for positive pions, kaons and protons in the 
most central 5 % of Au+Au collisions for y = at ^Jsnn = 130 GeV (51). 
Similarly, in Figure 5 STAR (47), PHENIX (53) and BRAHMS (58, 59) data 
are shown at ^/s^n = 200 GeV. Note the scaling factors 10 and 1000 for kaons 
and protons, respectively. An important issue of uncertainty in the calculation 
is the dependence of the results on the decoupling temperature. This is shown 
by plotting the results for freeze-out temperatures T^ ec = 150 MeV (solid lines) 
and 120 MeV (dotted lines). Note that the normalization of the pion spectrum 
is almost independent of the decoupling temperature. Because pions provide the 
main contribution to the total multiplicity, dN tot /dy depends only weakly on 
Tdec- However, the multiplicity of heavier particles is very sensitive to T^ ec , as 
expected. The shape of all three spectra changes clearly with the decoupling 
temperature. This follows from the increase in flow velocity, characterized by the 
increase of the effect with the increasing mass of the particle. The aim of this 
calculation has not been to find a best fit to the data, but the results show that 
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both the normalization and the slope of the data at momenta p T < 3 GeV can 
be described quite well with a single T^ ec in the neighborhood of 150 MeV. 
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Figure 4: Transverse momentum spectra 
of positive pions, positive kaons, and pro- 
tons at y — in the most central 5 % 
of Au+Au collisions at ^/snn = 130 GeV. 
The solid and dotted lines show our hydro- 
dynamic results for freeze-out temperatures 
T dcc = 150 MeV and T dcc = 120 MeV, re- 
spectively. The PHENIX data (51) is plot- 
ted with the given total error bars. Note 
the scaling factors 10 and 1000 for kaons 
and protons, respectively. Both the hydrody- 
namic result and the PHENIX data contain 
the feed-down contributions from hyperons. 
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Figure 5: As Figure 4 but at y/s NN = 
200 GeV. The PHENIX data (53) and 
the BRAHMS data (58, 59) are shown 
with statistical errors and the STAR 
data (47) with the given total error 
bars. The hydrodynamic calculation and 
the PHENIX data are without the hy- 
peron feed-down contributions whereas 
the STAR and BRAHMS data contain 
the feed-down. 



For the identified particles in Figures 4 and 5, the measured spectra do not 
extend to large enough pt to show clearly the deviation from the hydrodynamic 
results, with the exception of proton spectra. If the decoupling temperature is 150 
MeV to reproduce the normalization, the slope tends to be too steep. The proton 
yield from jet fragmentation, as explained in detail in Reference (44), does not 
seem to be large enough in the pr ~ 3 . . . 5 GeV region to bring the calculation 
into agreement with the data. In Figure 6, the spectra of two other heavy particles 
are shown, those of antilambdas and antiprotons. These show the same trend 
as protons, pointing to the need for separate chemical and kinetic decoupling 
when describing simultaneously the details of all spectra. Studies with separate 
chemical and kinetic decoupling, in which the stable particle numbers are fixed 
after chemical freeze-out, indeed show that the spectra of pions and kaons become 
almost independent of the kinetic decoupling temperature Id e c,kiin whereas the 
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(anti)proton spectra widen with decreasing Td ec ,km (40, 41). There are, however, 
claims in the literature that separate chemical and kinetic decoupling lead to a 
worse overall fit to the slopes of pt distributions than what can be achieved by 
requiring chemical equilibrium until kinetic freeze-out (60, 61). Studies exploring 
the effects of initial time, the shape of initial distributions, and the value of 
ec,chem while using two separate freeze-outs are needed to settle the issue. 
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Figure 6: Transverse momentum spectra 
of antiprotons and antilambdas at y = 
in the most central 5 % of Au+Au colli- 
sions at y/SNN — 130 GeV. Our hydro- 
dynamic results are for freeze-out temper- 
atures Tdcc = 150 MeV (solid line) and 
Tdoc = 120 MeV (dotted line), with hy- 
peron feed-down contributions included, 
as in the PHENIX p (51) and A (52) data 
and the STAR data (48). 
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Figure 7: Transverse momentum spec- 
tra of charged particles at r\ = (av- 
eraged over \r}\ < 0.1) in the most cen- 
tral 5 % of Au+Au collisions at ^snn = 
200 GeV. Our hydrodynamic results are 
shown for Tdcc = 150 MeV (solid line) 
and Tdcc = 120 MeV (dotted line). The 
pQCD fragmentation results are shown 
with (shaded band) and without (dashed 
line, see the text) energy losses. The data 
is taken by STAR (45), PHENIX (49), 
PHOBOS (55), and BRAHMS (57). 



The range where hydrodynamics can be used to describe the hadron spectra 
is indicated clearly in Figure 7, which shows results from our hydrodynamical 
calculations and from a pQCD jet calculation, followed by an energy loss in 
the medium before the jet fragments into hadrons (see below). The STAR and 
PHOBOS data are plotted with the given total error bars, the PHENIX data 
by adding the given statistical and systematic errors in quadrature, and the 
BRAHMS data with the given statistical error bars. 

The transverse spectrum up to ~ 3 GeV is similar to that of the dominant pion 
component shown separately in Figure 5. It has the shape typical of the spectrum 
from hydrodynamic calculations, falling off roughly exponentially. In the region 
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Pt ~ 3 . . . 4 GeV, there is a large change in the slope, indicating a change in the 
overall production mechanism. The calculation of primary production proceeds 
through the hard and semihard interactions between the partons of the incoming 
nuclei, both in the case of initial conditions for hydrodynamical equations and the 
energy loss of the jets. In the hydrodynamic calculation, the produced partons 
are assumed to thermalize quickly and then undergo hydrodynamic expansion 
in local thermal equilibrium until the freeze-out. In the calculation of jet frag- 
mentation after energy loss, the produced high-energy partons are assumed to 
survive the thermalization, but lose energy when traversing the thermal medium 
formed by the lower-energy partons. When energy loss and fragmentation are 
taken into account, the original energy of the partons that fragment to hadrons 
°f Pt 3 GeV must be of the order of ~ 6 GeV or greater. It turns out that 
the contribution of partons with pt > 4 GeV to the production of (transverse) 
energy is less than 5 %, justifying as a good approximation the assumption that 
all produced transverse energy is thermalized. The details of the jet-energy-loss 
and fragmentation calculation are explained in References (44) and (62). 
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Figure 8: Predictions for the transverse momentum spectra of charged hadrons at T) = 
in the most central 5 % of Pb+Pb collisions at Large Hadron Collider (LHC) energy 
\J snn — 5500 GeV. The shaded band in the hydrodynamic results shows the freeze- 
out temperature interval 120 . . . 150 MeV. The solid curve labeled "pQCD" corresponds 
to the pQCD fragmentation results without energy losses. The shaded band labeled 
"pQCD + E-loss" describes the uncertainty in the pQCD fragmentation results with 
energy losses. 

Although the two-component approach, hydrodynamically expanding thermal 
matter as the source of low-pT hadrons and jet fragmentation after energy loss 
producing the high-pr hadrons, seems reasonable and justified, adding them 
straightforwardly would be too naive. In the region of turn-over from one mech- 
anism to the other, both contribute, but part of the hadrons cannot be assigned 
to either component. Other mechanisms like recombination can also contribute 
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in this region (63). Accurate data in this area would be useful in understanding 
both the energy loss and thermalization of produced partons. 

As an example of the dependence on the collision energy, ^/sjvtv, Figure 8 shows 
an extension of the calculation to the CERN LHC, with ^/snn = 5500 GeV for 
the heavy ion collisions. In the calculation of initial conditions from the primary 
parton interactions, all parameters are fixed except for the collision energy. The 
saturation scale changes from p sat = 1.16 GeV at RHIC energy to p SSL t = 2.03 
GeV at the LHC. The total multiplicity increases from dNtot/dy « 1000 at RHIC 
to 4500 at the LHC. The initial thermal densities are higher and lead to longer 
expansion and stronger transverse flow at the decoupling. This is seen in the 
change of the region where the component of thermal particles goes over to the 
component of particles from jet fragmentation with energy loss. At RHIC the 
transition is centered around px ~ 3 GeV, whereas at the LHC it is predicted to 
be pt ~ 5 GeV. The larger pt region where thermal particles dominate should 
also be seen in elliptic flow. At the present RHIC energies, the hydrodynamic 
predictions of elliptic flow start to overshoot the data above pr ~ 2 GeV, but if 
the thermal component grows as predicted in Figure 8, the elliptic flow parameter 
V2 should follow the hydrodynamical calculation to larger transverse momenta. 



3.3 Elliptic Flow 

The particle production in primary collisions is azimuthally isotropic, whereas 
the reaction zone in noncentral collisions is not, but has an elongated shape. If 
produced particles rescatter, the particles moving in the direction of the longer 
axis of the reaction zone are more likely to change their direction than the parti- 
cles moving in the direction of the shorter axis. Therefore the observed emission 
pattern of particles will be azimuthally anisotropic, and the more frequent the 
rescattering, the more anisotropic the particle distribution. 

In this way, the anisotropy of the final particle distribution is a measurement 
of the frequency of rescatterings during the dense phase of the collision. This 
anisotropy can be quantified as the coefficients of the Fourier expansion of the 
azimuthal particle distribution (64): 

dN dN 

(1 + 2vi cos(0 - 4> R ) + 2v 2 cos 2(0 - <j> R ) + • • • ), 



dyd(j), p 2irdy 
dN dN 
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where (j) R is the azimuthal angle of the event plane (the plane spanned by the 
beam direction and the impact parameter). Assuming that the experimental 
uncertainties in event-plane reconstruction can be corrected for, each event can 
be rotated such that (f) R = 0. The first and second coefficient of the expansion, v\ 
and V2, are usually referred to as directed and elliptic flow, respectively. Because 
the system is usually thinner in the direction parallel to the impact parameter, the 
in-plane direction, than in the out-of-plane direction, the value of v 2 is positive. 
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At midrapidity, all uneven coefficients are zero owing to symmetry. At SPS 
and RHIC energies, the directed flow, v\, is expected to be very small and most 
of the experimental and theoretical interest has been directed toward measuring 
and analyzing the elliptic flow, v<i- Recently the higher harmonics, v^v§ and v$, 
have also been measured (65-67). 

In a hydrodynamic picture, the buildup of momentum anisotropies is easy to 
understand in terms of pressure gradients. The average pressure gradient between 
the center of the system and the surrounding vacuum is larger in the direction 
where the collision system is thinner. Therefore, the collective flow is stronger 
in that direction and more particles are emitted there than in the orthogonal 
direction, where the collision system is longer. 

As mentioned above, the more the particles rescatter, the larger the observed 
anisotropy. Because hydrodynamics assumes practically infinite scattering rate 
and zero mean free path, it is often assumed to give an upper limit of anisotropy 
at fixed impact parameter (25). However, this upper limit depends on the EoS 
and, in principle, it is possible that hydrodynamical description with very soft 
EoS would give a smaller anisotropy than, for example, a microscopic cascade 
description (68). 

If the freeze-out happens at the same temperature for all particle species, a 
signature of hydrodynamic flow is that the heavier the particle, the flatter the 
slope of its pT-spectrum. Similarly, the py-averaged elliptic flow vi increases when 
particle mass increases. However, the pT-differential elliptic flow V2(pr) has the 
opposite behavior: The heavier the particle, the smaller the anisotropy at fixed 
Pt- The apparent discrepancy has a simple explanation: v<i is not an additive 
quantity, but when pr-averaged ^2 is calculated from p^-differential V2(pr), the 
latter is weighted by the particle distribution. Thus, the flatter px distribution 
of a heavier particle weights more the high-pr region, where V2(pr) is larger. 
Therefore, even if V2{pr) is smaller at all px for a heavier particle, the pr-averaged 
V2 can be larger than V2 for a light particle. Whether this happens in practice 
and how large the differences are depend on the details of the flow profile, i.e. 
expansion dynamics, and the resonance decays. 

The mass ordering of the low-py anisotropy has its origin in the behavior of 
boosted particle distributions. As is well-known, transverse flow shifts the pr- 
distributions to larger values of pt- In the extreme case in which the speed of the 
collective motion is the same everywhere, as in the case of a thin shell expanding 
with a velocity v, the particle distribution develops a maximum at some finite 
Pt (the so-called blast wave peak (69)) and a local minimum at pr = 0. The 
heavier the particle, the larger the pt where the distribution peaks. Compared 
with a case without transverse flow, the particle yield is thus depleted at low 
Pt- The heavier the particle and the larger the flow velocity vt, the larger the 
depletion. Correspondingly, if the flow velocity is larger in the in-plane than in 
the out-of-plane direction, the low-p? depletion is larger for particles moving in 
the in-plane direction than the out-of-plane direction and the overall anisotropy, 
V2, is reduced. This reduction and the range in which it occurs increases with 
the particle mass and average transverse flow velocity. In the extreme case of 
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Figure 9: Elliptic flow V2 of charged hadrons as a function of centrality at (left panel) 
Pb+Au collisions at ^/snn — 17.3 GeV measured by CERES (71) and at (right panel) 
Au+Au collisions at ^snn = 200 GeV measured by STAR (66). Hydrodynamical 
calculations are taken from References (31, 71). Note that the hydrodynamical results 
are not directly comparable because of different pt cuts implemented. 



a thin expanding shell, this reduction can be so strong that it reverses the sign 
of the anisotropy and t>2 becomes negative. When the thin shell is replaced by 
a more realistic velocity profile, the peak in transverse-momentum distribution 
disappears. Similarly, a more realistic velocity distribution weakens the reduction 
of V2 at low pt, but the mass ordering of v-i at low pt remains. Whether some 
particles depict positive or negative V2 at low px depends on the details of the 
flow velocity of the source. 

For relativistic pt > m, the particle mass does not play any role in the thermal 
distribution and, consequently, V2(pr) of different particles converge. In a simple 
model in which the flow- velocity profile is replaced by its average value, V2(pt) 
increases with pt and approaches unity asymptotically. The details of the flow- 
velocity profile can change this behavior, but so far no hydrodynamical calculation 
has reproduced the experimentally observed saturation of elliptic flow. 

3.3.1 CENTRALITY DEPENDENCE 

In a hydrodynamic description, the final anisotropy of particles is almost directly 
proportional to the geometrical anisotropy of the initial state (70). The propor- 
tionality is, however, nontrivial and depends on the applied decoupling criterion. 
When the impact parameter increases and the collision becomes more periph- 
eral, the collision system becomes more and more elongated and its geometric 
anisotropy increases. We can thus expect the observed momentum anisotropy to 
increase as well. The data in Figure 9 shows increasing elliptic flow with decreas- 
ing centrality both at SPS and RHIC energies, but the magnitude of the flow 
differs from the hydrodynamical result. At SPS energy the data is consistently 
below the calculation, whereas at RHIC energy {^Jsnn = 200 GeV in Figure 9) 
the data is reproduced up to semi-central collisions but is below the calculation 
at peripheral collisions. 
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The failure of hydrodynamics to describe the anisotropy in most peripheral 
collisions and at SPS energy is often explained by a lack of necessary thermaliza- 
tion owing to the small size and/or particle number of the collision system (72). 
An alternative explanation assumes that the initial partonic state is sufficiently 
thermalized, but the final hadronic state has such a large viscosity that it cannot 
be modeled using ideal-fluid dynamics (73). The latter approach has been tested 
using so-called hybrid models in which the plasma phase and phase transition are 
described using ideal hydrodynamics, but the hadron phase is described using a 
cascade model (61, 73-75). The centrality dependence of py-averaged elliptic 
flow at RHIC has been reproduced nicely using such a hybrid approach (61, 73), 
but the results for px differential vi at ^/snn = 200 GeV collision energy are not 
available at the time of this writing. 

At the beginning of this section it was argued that hydrodynamics leads to the 
largest possible anisotropy. However, at most central collisions at RHIC, the data 
tends to be above the hydro dynamical calculation. Fluctuations in the shape of 
the initial system may explain this. Owing to these fluctuations, the initial shape 
of some events in almost-central collisions can be in-plane elongated, even if the 
shape on most events is out-of-plane elongated. Thus elliptic flow is negative in 
some events, but because experimental analysis measures the magnitude of the 
anisotropy, not its sign, elliptic flow is measured as positive in all events and the 
measured value is larger than the average value. The initial state of hydrody- 
namical calculation, however, is an average initial state in which fluctuations of 
the spatial anisotropy cancel each other and the calculated anisotropy is smaller 
than measured (76). Preliminary calculations in which the initial-state fluctua- 
tions are included favor this interpretation by leading to better reproduction of 
the data (77). 

The general trend is that a stiffer EoS and a lower freeze-out temperature lead 
to larger pr-averaged how if nothing else in the model is changed. This also 
changes the single particle distributions. If these are still required to fit the data, 
additional changes are required. For example, a stiffer EoS usually necessitates 
a higher freeze-out temperature. The combined effect largely cancels and the 
final pT-averaged anisotropy is almost unchanged in semi-peripheral collisions in 
which a hydrodynamical description works best (31). 

3.3.2 MOMENTUM and PARTICLE SPECIES DEPENDENCE 

Hydrodynamic calculations at SPS and RHIC energies lead to anisotropy, which 
increases with increasing px and approaches unity asymptotically. Simple parame- 
trizations of flow (so-called blast-wave models) also lead to this kind of behavior, 
which differs from experimental observations in which V2 saturates at high px- 
At mid rapidity, the agreement between the data and calculations depends on 
energy in the same way as for the centrality dependence: At SPS energies the 
calculations overestimate the data, whereas at RHIC energy a good agreement 
can be reached. The px range in which the data can be reproduced depends on 
the particle species. Charged hadrons and pions can be fit up to px ~ 1.5 GeV, 
whereas protons follow the calculations up to px = 2.5-3 GeV in minimum-bias 
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Figure 10: Pion (left panel) and antiproton (right panel) «2(pt) at midrapidity in 
minimum-bias Au+Au collisions at ^/snn — 200 GeV. Hydrodynamical results are 
labeled CE for chemical-equilibrium result (31), PCE for chcmical-nonequilibrium re- 
sult ((40, 61); T. Hirano, personal communication) and H+C for hydro+cascade hy- 
brid model at ^/s^n = 130 GeV energy (73). The data are taken by STAR (78) and 
PHENIX (79). 



collisions. This behavior is qualitatively similar to that seen in Figures 4 and 5 
for transverse momentum spectra, where hydro dynamically calculated spectra fit 
the data up to pt ~ 3 GeV. 

A hydrodynamic description predicts a characteristic mass dependence of el- 
liptic flow at low pt- The higher the mass, the lower the V2- How large this 
difference is depends on the details of the flow profile and therefore on the EoS. If 
chemical equilibrium in the hadronic phase is assumed, the differential anisotropy 
of pions can be well reproduced when the pt spectra of pions is reproduced. In 
such a case, the proton V2(pt) depicts sensitivity to the phase transition. If the 
phase transition takes place in a narrow temperature interval and has large latent 
heat, the proton differential anisotropy can almost be reproduced. If there is no 
phase transition, the calculated proton anisotropy is clearly above the data (31). 
The status of strange particles is less satisfactory. Kaons and A-baryons show 
similar dependence on the EoS than protons, but the difference between the data 
and hydrodynamical calculation is larger than in the case of protons (31). 

Unfortunately, it is not yet possible to use the apparent sensitivity of proton 
V2{pt) to the EoS to quantitatively constrain the EoS. If the requirement of 
chemical equilibrium is relaxed and one uses separate chemical and kinetic freeze- 
out temperatures, the fit to pion i>2(pr) is lost (see Figure 10 and Reference (40)). 
However, if the hadronic phase is described using the RQMD transport model, 
as in Reference (73) for the ^/snn = 130 GeV collision energy, the yields are 
correct and the V2{pt) is described as well as in the case of chemical equilibrium. 
Thus, there is a considerable uncertainty in the description of the hadronic stage 
of the evolution, which makes it impossible to draw final conclusions about the 
EoS needed to describe the differential anisotropy. 

Another uncertainty here is the effect of viscosity. The ability of ideal-fluid hy- 
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drodynamics to reproduce the vi data at RHIC has been interpreted to mean that 
the value of shear viscosity in QGP is particularly low (60, 80). However, viscos- 
ity has been estimated to decrease elliptic flow (4) , and chemical non-equilibrium 
increases elliptic flow (40, 41). Hirano h Gyulassy have argued that the plasma 
has sufficiently low viscosity to allow its modeling using ideal hydrodynamics, 
but the dissipative effects in the hadronic phase are the reason for the failure of 
chemically non-equilibrium hydrodynamics to describe the data (60). This claim 
seems to be validated by the ability of hybrid models to describe the data, but 
final conclusions must wait for the complete results of hybrid calculations and 
the results of viscous calculations (5, 81). 

3.3.3 PSEUDORAPIDITY DEPENDENCE 

The pT-averaged elliptic flow has quite a strong dependence on pseudorapid- 
ity (78, 82-84). In a narrow region around midrapidity, \r\\ < 1, elliptic flow re- 
mains approximately constant (78, 85) but decreases strongly towards larger pseu- 
dorapidities. However, charged particle multiplicity depicts much wider plateau 
around midrapidity than elliptic flow (86). 

The purely hydrodynamical calculations have not reproduced the pseudora- 
pidity dependence of elliptic flow, although the same calculations reproduce the 
multiplicity as a function of pseudorapidity (40, 87). However, a hybrid model 
that reproduces the centrality dependence of elliptic flow also gives a reason- 
able description of its pseudorapidity dependence (61). The failure of ideal-fluid 
hydrodynamics has been interpreted in the same way as in the case of central- 
ity dependence - either as an incomplete thermalization of the system at large 
rapidities from the beginning (72) or as an effect of viscosity and incomplete 
thermalization at the late stage (61). There are also other open questions in the 
hydrodynamic treatment that may affect the results: The initial shape of the 
system, large deviations from boost-invariant flow, and different thermalization 
time and freeze-out temperature at different rapidities could all affect the final 
anisotropy and are mostly unexplored. It is therefore possible that the thermal- 
ized region where hydrodynamics works at RHIC energy is relatively narrow in 
rapidity, but final conclusions cannot be drawn yet. 

3.4 Two- Particle Bose- Einstein Correlations 

Information about the space-time structure of the system formed in a heavy-ion 
collision can be obtained by measuring the low-momentum correlations of identi- 
cal particles. For bosons these correlations are called Bose-Einstein correlations 
and the method for their interpretation is called HBT interferometry according 
to the originators of this method (88). Here we show only the basics of the HBT 
formalism as applied to hydrodynamical models and its most important results. 
A detailed explanation about this technique can be found in Reference (89), and 
the present status is discussed in recent reviews (90, 91). HBT in hydrodynamical 
context is more throughly discussed in Reference (92). 

Intensity interferometry is based on an analysis of the two-particle momentum 
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correlation function, 



^i^2 d ap id a p2 



g(g,g)= d ;^% (ii) 

that is, the ratio of a two-particle distribution and a product of two one-particle 
distributions. The correlator is usually written in terms of the momentum dif- 
ference between the two particles, q = Pi — P2, an d their average momentum, 
K = \{vi+Pi)- If the particles are emitted independently ("chaotic source") 
and propagate freely from the source to the detector, the two-particle distribu- 
tion is not equal to the product of one-particle distributions. At small values of 
relative momentum q, it is larger than the product of one-particle distributions 
owing to quantum statistical (wave-function symmetrization) effects. 

If there are no final-state interactions (or the spectra are corrected for them), 
the two-particle correlator C(q,K) is related to the emission function S(x,K): 

2 



C(q,K) w 1 + 



Jd 4 xS(x,K)e iqx 



(12) 



fd 4 xS(x,K) 

The emission function S(x, K) is the Wigner phase-space density of the emit- 
ting source. In the derivation of Equation 12 the emission function is assumed 
to be sufficiently smooth, i.e. S(x,K) « S(x,K + ^q) (see Reference (89)). 
Because both p\ and P2 are on-shell, the average momentum K is, strictly 
speaking, off-shell. In practice, however, on-shell approximation for K is used: 
K « ^/k 2 



m 



2 



In the hydrodynamic approach the quantum-mechanical Wigner phase-space 
density is replaced by a classical phase-space density at the time of freeze-out. 
When Cooper-Frye formalism is applied, it is given by 

Six K)~ 9 f da,(x')K^(x - x>) 

K ' ] (2tt)3 J exp{[K-u(x')-fi(x')]/T(x')}±l [ ' 

It is not possible to define uniquely the source function S(x, K) from the mea- 
sured correlation function C(q,K). The experimental data of two-particle cor- 
relations are therefore presented using some ansatz for the the source function. 
Usually this is done using a Gaussian form for the correlator. If the collision 
system is further approximated to be boost-invariant, the correlator for central 
collisions can be written in a particularly simple form in terms of three HBT 
radii: 

C(q, K) « 1 + exp[-R 2 (K T )q 2 - R 2 s {K T )q 2 s - Rj(K T )q 2 ]. (14) 

In this so-called Bertsch-Pratt parametrization, the coordinate directions are de- 
fined in such a way that out- (R Q ) and long-direction (R[) are parallel to Kt 
and beam, respectively, whereas the side-direction (R s ) is perpendicular to both 
Kt and beam. In a boost-invariant approximation, the radii depend only on the 
magnitude Kt because the particle emission in central collisions does not depend 
on the azimuthal angle (p and boost-invariance means that there cannot be any 
rapidity dependence. 

These radii do not correspond to the actual physical size of the source, but 
rather characterize so-called regions of homogeneity, the regions where particles 
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Figure 11: Pion source radii in Au+Au collisions at ^snn = 200 GeV. Hydrodynami- 
cal results are labeled CE for boost-invariant chemical-equilibrium result (26), PCE for 
nonboost-invariant result with separate chemical and kinetic freeze-outs ((40, 61); T. 
Hirano, personal communication) and H+C for hydro+cascade hybrid model (94). The 
data are taken by the STAR (95) and PHENIX (96) collaborations. 



with particular px are most likely emitted. For a Gaussian source, the HBT 
radii measure the following different combinations of space-time variances of the 
system (92): 

R 2 S (K T ) = (x 2 s ){K T ) (15) 
R 2 {K T ) = ((x -(3 ± t) 2 }(K T ) (16) 
R 2 {K T ) = {x 2 )(K T ), (17) 

where = Kt/K° is the transverse pair velocity, and space-time coordi- 
nates x are defined as distances from the "effective source center" x^{Kt) = 
x^ — (x^)(Kt), where brackets denote weighted averages over the source function 

The radii are thus independent of the actual coordinates of the emission, but are 
sensitive to variances of the geometry. 

Hydrodynamic calculations for RHIC energies predicted that a phase transition 
from a plasma to a hadron gas would increase the lifetime of the system (93). 
The long lifetime would increase the (3±t term in Equation 16 and thus increase 
the -R -radius and make the ratio R Q /R S large. The experimental data, however, 
shows no sign of this kind of effect and yields a ratio of R Q /R S w 1. 

Figure 11 shows some of the hydrodynamic calculations for the HBT radii at 
RHIC energy. Hydrodynamic calculations with assumed boost-invariance, chem- 
ical equilibrium, and a first-order phase transition (CE, solid line) tend to lead 
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to a too small sideward radius R s , too large outward and longitudinal radii R Q 
and Ri, respectively, and especially to a too large ratio R /R s (26). Both R Q and 
Ri can be made smaller if the system decouples sooner, i.e. in higher tempera- 
ture, but such an approach leaves R s basically unchanged and distorts the single 
particle spectra (26, 92). To a lesser extent, the same effect can be achieved by 
decreasing the initial time or increasing the transverse flow by a non-zero initial 
velocity field, but neither of these approaches changes the too small R s (26). 

Another way to reduce the longitudinal radius Ri is to relax the boost-invariant 
approximation (40, 97). When this approach is used with an EoS with separate 
chemical and kinetic freeze-outs (PCE, dashed line), i?; is close to the data. Both 
R and R s move closer to the data but are still too large and small, respectively. 
This approach also leads to problems with the elliptic anisotropy (see chapter 3.3). 

An approach that brings R s close to the data chooses a wide but flat initial 
distribution (97), which leads to slower buildup of flow from an initially larger 
source. In that case, the problem is again R Q , which is too large. This is expected 
because R Q is sensitive to the lifetime of the system, which becomes relatively 
large in this approach. 

One way to reduce the lifetime of the system and thus R Q is to change the EoS. 
As mentioned, one of the suggested signatures of a first-order phase transition is a 
long lifetime and large R Q . If one uses an EoS with a smooth crossover instead of 
a first-order phase transition, R Q decreases and R s increases (98). Unfortunately, 
even in that case, no good fit to the data is achieved. It is also questionable how 
this change in EoS would affect the elliptic anisotropy (31). 

Grassi et al. have suggested that the discrepancy between the data and calcula- 
tions is due to too simple a treatment of freeze-out on a sharp hypersurface, and 
a more realistic continuous emission of particles would lead to better results (99). 
However, when this is accounted for effectively in hybrid models in which the 
hadronic stage is described using a cascade transport model, the results are even 
worse (94). The particles are emitted from larger, longer-lasting volume than 
in a simple hydrodynamic description, and correspondingly, R s is larger and re- 
produces the data (H+C, dotted line in Figure 11). Unfortunately, the longer 
lifetime also leads to an even larger R Q . 

Another possible reason for the discrepancy between the data and calculations 
is viscosity (4, 100). Initial calculations (5) show that it has the desired effects, 
but whether they are large enough remains to be seen. Again, the effect of 
viscosity on elliptic flow is large, and it is unknown if a viscous model could 
reproduce both the HBT radii and elliptic anisotropy. 

4 ELECTROMAGNETIC EMISSION 

All the observables described in the previous section 3 are hadronic observables. 
By definition, the hadrons of the system interact with each other, and the dis- 
tributions and yields of hadrons are fixed late in the evolution of the system, 
when interactions cease, the distributions and yields freeze out, and the particles 
decouple. Therefore, those observables characterize the properties of the particle- 
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emitting source at the end of the system evolution, but not the history of the 
system during the evolution. In principle, it is possible to have very different 
dynamics producing similar final states. 

Possible observables that are sensitive to the entire evolution of the collision 
system are photon and lepton-pair distributions. Because these particles interact 
only electromagnetically, their mean free paths are much longer than those of 
hadrons. They can thus escape the collision system without rescattering and 
carry information about the conditions in which they were formed. However, the 
photon and dilepton spectra get contributions from all stages of the evolution, 
which makes it difficult to disentangle the signal coming from the hot, dense 
stage of the collision. To describe the different contributions to electromagnetic 
spectra, we follow the terminology of Reference (101): 

• prompt photons and leptons are produced in the primary collisions of in- 
coming partons; 

• thermal photons and leptons are emitted in the collisions of quarks and 
gluons during the plasma phase and in the collisions of hadrons in the 
hadronic phase; 

• decay photons and leptons are decay products of hadrons; 

• direct photons and leptons are the sum of prompt and thermal photons. 

The thermal photon production depends strongly on temperature via the factor 
exp(— pt/T). Therefore the early stage, when the matter is hottest, should dom- 
inate the photon emission, and the measurement of photon spectra should be 
an effective thermometer for the temperature achieved in the collision. However, 
prompt photons follow a power-law distribution p^ n and dominate at high pr- 

The hydrodynamic model can be used to calculate the thermal and decay 
contributions, but the prompt photons and leptons require a separate pQCD 
calculation. The calculation of decay photons is relatively straightforward, and 
it proceeds in the same way as the calculation of hadron spectra from resonance 
decays described in section 3.1. One calculates the distribution of hadrons at 
freeze-out and applies the relevant decay kinematics and branching ratios to get 
the spectra of decay photons and leptons. To calculate the thermal yield, the 
production rate of photons or leptons in a thermal system, dR/d 3 p(E, T, fx), has 
to be integrated over the space-time volume of the system: 

„dN f l4 r , , dR QGP , m , 

~dJp = J I w{e,p B ) E—p—(p-u,T,Li B ) 

dR UG ^ 
+ [l-w(e,pB)]E-p-(p-u,T,ti B )l , (19) 

where the production in a plasma and in a hadron gas (HG) is written separately. 
The factor w(e,pb), which expresses the volume fraction of plasma, is unity in 
the plasma phase, zero in HG and between unity and zero in a mixed phase. 
Hydrodynamics provides the space-time evolution of the system, whereas the 
production rates in thermal matter are an input independent of hydrodynamics. 
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In the plasma phase, the photon production is dominated by the QCD Compton 
and annihilation reactions, qg — > (77, qg — > 57, and — > 37. In lowest order, the 
production rate due to these processes was calculated in References (102, 103). 
However, some formally higher-order processes are strongly enhanced by collinear 
singularities and also contribute to order a s (104, 105). The resummation of these 
contributions was shown to be possible and was carried out by Arnold et al. (106), 
completing the order a s analysis of the photon-emission rate. A parametrization 
of the rate was also given in Reference (107). 

The calculation of the photon-production rate in a hot hadron gas is less com- 
plete than the rate in a plasma owing to the multitude of different hadron species 
and photon-producing interactions and owing to the model dependence of the 
calculations (see References (101, 108)). The standard rate in the literature is 
the one calculated in Reference (102), where photon production in scattering and 
decay processes nir — ► pj, np — > n^f, to — > 7r7, and p — ► tttt^ was calculated using 
pseudoscalar-vector Lagrangian with coupling constants determined from free p 
and uj decays. These rates are often supplemented with a production rate via 
ai-mesons from Reference (109). 

The role of different channels in photon production was further studied using 
chiral Lagrangians (110). Unfortunately, it was not possible to fix the model 
parameters unambiguously in this work, which led to a factor of three uncertainty 
in the final rates. In the context of dilepton production, it was later possible to 
fix the model parameters much better (111), and this approach was used in a 
recent calculation by Turbide et al. (112). In that work, the study was extended 
to cover photon emission from heavier meson resonances, strange particles, and 
baryons. Another recent analysis of photon production from a hadron gas was 
done by Haglin (113), who studied the effect of strange particles and higher-order 
processes achieving a rate larger than the standard rate (102) by a factor of two 
at large qx and by an order of magnitude at low qx- 

There are still uncertainties in the calculation of the photon-production rate in 
hot hadron gas. Surprisingly, even after all the improvements in the calculations, 
the statement made in Reference (102) is still valid. At the same temperature, 
the production rate per unit volume in a plasma and a hot hadron gas is approx- 
imately equal, and they both "shine as brightly" . However, the emission rate per 
unit entropy is larger in hadron gas. 

The main contribution to dilepton production in plasma comes from the an- 
nihilation process qq — ► 11. The rate calculated in lowest order in a baryon-free 
plasma can be found in textbooks (114) and was calculated for finite baryon 
chemical potentials in Reference (115). At small values of invariant lepton-pair 
mass, corrections of the order aa s to this rate become important (116), but in 
heavy ion collisions lepton pairs from Dalitz decays of final mesons produce a 
larger background (117). Multi-loop calculations similar to those done to cal- 
culate the photon rate have also been carried out for high-p^ pairs with small 
invariant mass (118). These calculations have resulted in rates somewhat larger 
than the first-order calculations. First attempts to calculate lepton production 
using a lattice-QCD formalism have also been done (119). The preliminary re- 
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suits are quite close to the perturbative rate, at least in some parts of the phase 
space. 

The observation of large excess dileptons in the mass region below the p- meson 
mass in Pb+Au collisions at = 17.3 GeV energy at the CERN-SPS (120) 

has fueled considerable theoretical interest in studying the lepton-pair emission 
in a hot hadronic gas (108, 121). The main problem of these studies has been 
whether and how the properties of mesons change in medium and how these 
changes are reflected in the rates. A rate calculated by Gale & Lichard (122) 
using free-particle properties is often used as a benchmark in comparisons with 
more sophisticated approaches. In the calculations of Rapp et al. (123) and 
Eletsky et al. (124), the basic assumption is that the spectral density of /5-meson 
changes in medium. These calculations are technically very different, but produce 
qualitatively similar rates (125). An alternative approach pursued by Brown & 
Rho (126) assumes that the p-meson mass decreases in the medium. 

4-1 Photons at SPS 

Direct photon production in y / sjv r /v = 17.3 GeV Pb+Pb collisions at the CERN- 
SPS was measured by the WA98 collaboration (127). Several authors have com- 
pared this data with hydrodynamical calculations (128-132). All authors agreed 
that the photon spectrum could be explained if one assumes sufficiently hot 
(T > 200 MeV) initial state, but the required initial temperature varied largely 
from T ~ 200 MeV (129) to T = 335 MeV (128). The large difference is owing 
mainly to different assumptions in the calculations. 

One factor that explains the largely varying initial temperature is the use of 
different rates in a hadron gas. Alam et al. (129) assumed in-medium modi- 
fications to hadron properties both in the EoS and in production rates, which 
enhance the photon emission at lower temperatures, allowing cooler initial state. 
The full order as rate for photon production in plasma was used only in the 
most recent paper (132), but at SPS energy the different rates in plasma cause 
significant differences in the final yield only at relatively large values oipT- 

The initial state of the system was also chosen in different ways in different 
calculations. Especially, the assumption of finite transverse flow velocity at the 
beginning of the hydrodynamic evolution leads to lower temperatures. Because 
the rates are proportional to exp(— p • u/T), where p is the four momentum of 
the photon and u is the flow four-velocity, stronger transverse flow allows lower 
temperatures to produce equal yield at high px- Peressounko & Pokrovsky (130) 
argued the necessity of such an initial flow, and Alam et al. (129) and Chaud- 
huri (131) later studied its effects. It can be argued that gradients in initial 
particle production would lead to buildup of flow during thermalization, but it is 
very difficult to quantify how large flow velocities could build up this way. Peres- 
sounko & Pokrovsky also argued that the pion spectra especially necessitates the 
initial flow, but the authors of this review have not been able to fit the hadron 
spectra if initial transverse flow is assumed. 

The hadron spectra were reproduced in References (130, 132), and, when no 
initial transverse flow is assumed, also in Reference (129). In the other two cal- 
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culations (128, 131), the initial state was only required to have the same entropy 
as the final-state particles. It is thus unknown whether these calculations are 
consistent with the hadron data. 

With the exception of Reference (132), boost-invariant hydrodynamics was 
used in these calculations. If high initial temperature is required, this assump- 
tion leads to short initial, i.e. thermalization, time, r ~ 0.2 fm/c (128). At 
SPS energy it can be argued that such a short initial time is ambiguous be- 
cause the longitudinal extension of the colliding nuclei is larger than ctq. This 
makes the application of boost-invariant expansion uncertain for times r < 1 
fm/c. In Reference (132) this problem was solved by using nonboost-invariant 
hydrodynamics, where longitudinal geometry is explicit and the initial time does 
not appear. Also, the ambiguity in choosing the initial state was studied and it 
was shown that several EoSs and initial states reproduced both the hadron and 
photon data (132, 133). 

To characterize the results, Figure 12 shows the photon spectrum calculated in 
Reference (132) and compares it with WA98 data. The calculation was done using 
two different equations of state (EoS A and H) and two different initial states 
(IS 1 and IS 2). EoS A contains a phase transition from hadron gas to quark- 
gluon plasma at T c = 165 MeV whereas EoS H is a purely hadronic equation 
of state. IS 1 has a very peaked initial density distribution in the longitudinal 
direction, whereas IS 2 has a flatter distribution (see Reference (133)) and smaller 
maximum temperature, which is more consistent with the assumption of hadronic 
EoS. In both cases, the hadron spectra are reproduced and, as shown in the figure, 
the calculated photon spectrum is within the experimental error bars. Thus, the 
conclusion of the hydrodynamic studies of photon emission at the SPS is that high 
temperature initial state is needed to reproduce the measured photon spectra, 
but a phase transition to plasma is not necessarily required. 

4.2 Photons at RHIC 

At the time of this writing, the situation of the photon data at RHIC is becoming 
very interesting: Preliminary data presented (135) at the latest Quark Matter 
2005 meeting indicates a clear excess of photons over decay and prompt photons 
in the transverse-momentum range up to ~ 3 GeV. Earlier measurements have 
been inconclusive owing to the large error bars, but the new method to extract 
the photon yield from the measurements of low-mass e + e~ pairs appears more 
promising, and if the preliminary results are confirmed by the full analysis, this 
data offers a long-sought direct probe into the earliest moments of the collision. 

Many authors have predicted the photon emission at RHIC and the LHC (101, 
105, 129, 136-141). Owing to uncertainties in the initial state, these predictions 
serve mostly as order-of-magnitude estimates, but they also address the question 
of whether the thermal photon yield would be larger than the prompt photon 
yield at any value of pr- The initial conditions of more recent studies have been 
constrained to produce the total hadron multiplicity, and in Reference (101), 
the thermal photons are compared with the calculated yields of decay photons 
both from thermal pions and prompt pions from jet fragmentation. The measured 
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Figure 12: The sum of thermal pho- 
ton emission at the SPS from hydrody- 
namic calculation (132) and prompt pho- 
ton emission (134) compared with WA98 
data (127). 



Figure 13: Thermal photon emission at 
RHIC from hydrodynamic calculation us- 
ing two different initial times. Solid lines 
labeled TOT indicate the total yield of 
thermal photons whereas dashed lines in- 
dicate emission from plasma. The pQCD 
calculation for prompt photons (101) and 
preliminary PHENIX data (135) are also 
shown. 



spectrum of tt s is compared with the calculations and the data is well described. 
The hydro dynamical calculations with the same initial conditions are also com- 
pared with other hadron data and the overall agreement is very good (44). The 
perturbative QCD calculation to NLO in all quantities entering the calculation 
is supplemented with the energy loss of produced jets in the thermal matter. 

The main conclusions from the studies in (101), supported for the hadron ob- 
servables by (44), are as follows: The understanding of the hadron spectra in 
terms of hydrodynamic and pQCD calculations is quite good. This means that 
the photons from hadron decays are well under control when comparing different 
sources in the calculations. The other main sources of photons are prompt pho- 
tons from primary interactions, including the photons from jet fragmentation, 
and the photons from secondary interaction in produced matter, the thermal 
photons 3 . These photon sources have quite distinct transverse-momentum de- 
pendence with the cross-over from thermal photons to prompt photons taking 
place at around pt ~ 3 GeV, the region where the behavior of the preliminary 
data also changes. 

The simplest hydrodynamic calculations assume a scaling expansion in the 
longitudinal direction and ignore the transverse expansion. The qfirst assumption 
can be argued to be reasonable in the central rapidity region because at RHIC 

3 More generally the photons from secondary interactions among the produced particles can 
originate also from non-thermal processes like a high-energy quark producing photons when 
Compton scattering from a lower energy thermal gluon (142). 
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energy the Lorentz gamma factor is ~ 100, indicating a time interval on the 
order of 0.1 fm/c for the nuclei to pass through one another. This is shorter than 
the shortest initial times used in the calculations. In the central rapidity region, 
the longitudinal components are not large and the acceleration of the longitudinal 
expansion is small, having little effect on the multiplicity density or the freeze-out 
time at y (143). Ignoring the transverse expansion cannot be justified, except 
for the photons emitted at the earliest times from the QGP. High-pr photons from 
QGP are almost insensitive to flow because they are emitted when the system 
is hottest (141), but the strong flow at the late hadronic stages enhances the 
emission of high-py photons from hadron gas (129). 

Predictions of the relative size of photon contributions at different values of px 
vary somewhat. For example, Srivastava (140) predicts that at RHIC energy, the 
photons from the QGP dominate at small values ofpr, i.e., px < 1 GeV, and the 
thermal photons at high px come mainly from the hadronic phase. More recent 
calculations lead to a conclusion that at RHIC multiplicity the contribution from 
plasma dominates for pt > 3, and at smaller transverse momenta the contribu- 
tion from plasma and hadron gas are the same size, with the latter slightly larger 
at smallest momenta (136, 138). The rates used in the calculations is one reason 
for the difference: In Reference (140) older rates that do not include all order 
of as terms are utilized, whereas the newer calculations are based on the full 
order of as (107) results. Also, note that the dependence of the plasma contribu- 
tion depends strongly on the assumed thermalization time, To- When comparing 
different predictions the first detail to be checked is To, see Reference (101). 

Although the pQCD calculation of prompt photons is not entirely under control 
at small transverse momenta owing to the uncertainty in the photon fragmen- 
tation functions (see discussion in References (101, 144)), the photon yield from 
secondary collisions, i.e. the yield of thermal photons, decreases more steeply and 
becomes negligible for px ^ 4 GeV. In the calculations thermal and prompt pho- 
tons become comparable at around px ~ 3 GeV, and because of the difference in 
the slopes, the uncertainty in where the contributions cross is not large. At small 
momenta below 3 GeV, thermal photons dominate (136, 138), but it is not clear 
whether this contribution is so large that it can be isolated from the pion decay 
background. 

Different contributions are compared with the preliminary photon data (135) 
in Figure 13. Here, initial times To = 0.2 and 0.6 fm/c are used in the hydrody- 
namical calculation, which correspond to average initial temperatures (T) = 340 
and 250 MeV, respectively. In the prompt photon calculation no intrinsic kx is 
included because the same calculation for photon production in p + p collisions 
at RHIC describes the data well. As indicated in Figure 13, the uncertainty in 
initial time causes more than an order of magnitude change in the thermal photon 
yield at high px- 

In most calculations for thermal photons at RHIC and the LHC, chemically 
equilibrated matter is assumed. However, it can be expected that the initial state 
is gluon dominated and quarks are suppressed compared to their equilibrium 
yields (145). This would lead to smaller emission rates at a given temperature, 
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but the suppression of quarks means effectively smaller number of degrees of 
freedom and a larger temperature for a given entropy. Detailed calculations (101, 
146) have indicated that these two effects largely cancel each other and that the 
final thermal spectra are quite similar in both equilibrium and non-equilibrium 
scenarios. 

4-3 Dilepton calculations 

Like the photon measurements, the dilepton mass spectrum at the SPS collisions 
measured by the CERES (120, 147) and NA50 (148) collaborations has been 
compared to hydrodynamical calculations several times (35, 125, 129, 133, 149- 
152). The low-mass dilepton (Mjj < m^) yield measured by CERES is dominated 
by emission from the late hadronic phase (see, e.g. Reference (151)) and constrains 
only the properties of the hadronic stage of the evolution. All these calculations 
agree that if meson properties in vacuum are used (122), the thermal yield is 
not quite sufficient to explain the observed excess. Thus, the experimental data 
seems to require modifications in meson properties, but so far the mass resolution 
has not been good enough to differentiate between a change in mass and changes 
in spectral density. Recently, there has been new preliminary data with better 
mass resolution (153, 154), but the conclusions are still being debated (155). 

Kvasnikova et al. (152) addressed the intermediate-mass (m^ < Mu < mj/y) 
dilepton yield measured by NA50. In their calculation they found that the excess 
in lepton pairs in this mass region could be explained by thermal emission in 
the same way as in the low-mass region. Even if the intermediate-mass region 
is expected to provide a window for observing the emission from plasma (156), 
they found that a modest contribution from plasma, ~ 20%, was enough to fit 
the SPS data. 

So far the genuinely hydrodynamic calculations of dilepton emission at RHIC 
have been rare, and different parametrizations for the space-time evolution of the 
system have been used instead (123). Medium modifications to meson properties 
depend on total baryon density, ptot = Pb+Pb- Because ptot at RHIC is essentially 
the same as at the SPS, the low mass dilepton spectrum at RHIC should show 
similar excess as seen at the SPS (123). 

At intermediate masses, the thermal yield is expected to be dominated by 
emission from plasma (157). However, owing to a larger cc production than at 
the SPS, the intermediate-mass dilepton yield can be dominated by correlated 
charm decays, unless the c-quarks rescatter significantly in the medium or can 
be identified and subtracted. Isotropization of c-quark momentum distributions 
would soften the dilepton mass spectrum, leaving a mass range in which thermal 
emission dominates (157, 158). An additional source of lepton pairs at RHIC 
is the interaction of jets with the surrounding dense plasma. According to re- 
cent calculations (159), the jet-plasma interactions may dominate over thermal 
dilepton emission at intermediate masses. 

At the time of this writing there are no calculations in which all these contri- 
butions are taken into account and are folded with a realistic time-evolution of 
the collision system. It will be interesting to see how the future data will look 



36 



P. Huovinen and P.V. Ruuskanen 



and if a contribution from plasma is needed to explain the data. The PHENIX 
collaboration has recently shown the first preliminary data of low-mass dileptons 
at RHIC (160), but the present experimental uncertainties are too large to draw 
any conclusions. These experimental shortcomings are currently being addressed 
by a detector upgrade (161). 

5 CONCLUDING REMARKS 

Hydrodynamics provides a well-defined framework to study many experimentally 
accessible features of a heavy ion collision. Some parts of the collision, such as 
the primary particle production of final-state matter, lie outside hydrodynam- 
ics, and some features of hydrodynamics, such as the freeze-out of final hadrons, 
have grave uncertainties. Nevertheless, a hydrodynamic description has robust 
features such as the the conservation laws, which are strictly enforced. The main 
assumption of hydrodynamics, the occurrence of frequent collisions in the final 
state, can also convincingly be argued for from the large observed multiplicity. 
Hydrodynamics describes the effects of collisions among the constituents, in par- 
ticular the momentum transfer between adjacent regions, in terms of pressure that 
arises microscopically from momentum transfer in the collisions. This should be 
a good approximation as soon as the momentum distribution of constituents is 
approximately isotropic. 

Hydrodynamics describes well the broadening and its mass dependence of 
hadron spectra resulting from the increase of transverse collective motion (flow). 
The collective motion can also be seen in the elliptic flow. At low transverse 
momenta, the observed elliptic flow can be described using hydrodynamics. Es- 
pecially, the observed mass ordering is typical for a hydrodynamic description. 
These are partly genuine predictions of hydrodynamics because the amount of 
initial production is fixed from the total multiplicity. Once the hadronic ob- 
servables are under control, the largest remaining uncertainty concerns the time 
scales of primary production and of (approximate) thermalization. Electromag- 
netic emission is very sensitive on these time scales, and the preliminary results on 
photon emission at RHIC may be the first indication that emission from the early 
moments of the collision can be resolved. The emission of lepton pairs around 
and below the <^>-meson mass, offers both a stringent test of the hydrodynamic 
description of the hadron phase and a tool to study the effects of medium on the 
properties of vector mesons. The amount of information that can finally be ob- 
tained, depends a good deal on the progress in experimental measurements, but it 
is likely that hydrodynamics will remain an important tool for phenomenological 
studies for a long time. 
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